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ABSTRACT 

This paper presents the results of an U band survey with FORS1/VLT of a large area in the tr Ori star-forming region.We combine the 
U-band photometry with literature data to compute accretion luminosity and mass accretion rates from the U-band excess emission 
for all objects (187) detected by Spitzer in the FORS1 field and classified by Hernandez et al. (2007) as likely members of the cluster. 
The sample stars range in mass from ~ 0.06 to ~ 1.2 M G ; 72 of them show evidence of disks and we measure mass accretion rates 
M acc between < 10~" and few 10~ 9 M Q /y, using the colors of the diskless stars as photospheric templates. Our results confirm the 
dependence of M acc on the mass of the central object, which is stronger for low-mass stars and flattens out for masses larger than 
~ 0.3 M Q ; the spread of M acc for any value of the stellar mass is ~2 orders of magnitude. We discuss the implications of these results 
in the context of disk evolution models. Finally, we analyze the relation between M acc and the excess emission in the Spitzer bands, 
and find that at M acc ~ 10~ 10 M Q /y the inner disks change from optically thin to optically thick. 

Key words. Stars: formation - Accretion, accretion disks - open clusters and associations: individual: <x Orionis 



1. Introduction 

In recent years, our knowledge of the properties of young stars 
in several star-forming regions has made enormous progress. In 
particular, Spitzer observations have provided new information 
on the IR properties of circumstellar disks, and the discussion on 
how disks evolve in time has gained new momentum. In addition 
to the classical viscous evolution, which dissipates disks by ac- 
creting their matter onto the central stars, other processes such 
as gravitational instabilities, photoevaporation by X-ray and UV 
radiation of the central star (Hollenbach et al. 2000; Gorti & 
Hollenbach 2009), and planet formation (Dullemond et al. 2007) 
have been recognized to be important. Photoevaporation and 
planet formation may both shape the SED (spectral energy dis- 
tribution) of so-called transitional disks, which have very low 
emission in the near- and mid-IR and strong excess emission 
at longer wavelengths (Calvet et al. 2002, 2005; D'alessio et 
al. 2006; Currie et al. 2009). In all cases, disks seem to come 
in all variesties: it is clear that neither time nor the mass of 
the system control them uniquely, and it seems very likely that 
disk properties and evolution depend also on the initial condi- 
tions, i.e., on the properties of the molecular core from which the 
star+disk system forms (Hartmann et al. 1998, 2006; Dullemond 
et al. 2006, Clarke 2007, Vorobyov & Basu 2009). 

An important contribution to this discussion comes from 
measurements of the mass accretion rate for well characterized 
samples of stars. However, there are only few systematic deter- 
minations of the mass-accretion rates in large samples of objects 
within the same star-forming regions, covering a large range of 
central masses, with well measured SEDs and complete to in- 
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elude also diskless stars, limited so far to p-Oph (Natta et al. 
2006) and Tr 37 (Sicilia-Aguilar et al. 2010). 

In this paper, we add a third region, <x Ori, to the list. The 
<x Ori cluster is ideally suited for this kind of study. It contains 
more than 300 young stars, ranging in mass from the bright, 
massive multiple system <x Ori itself (the spectral type of the 
brightest star is 09.5V, Caballero. 2007) to brown dwarfs. It is 
located at a distance of ~ 360 pc (Hipparcos distance 352^ 6 pc 
for the 09.5V star Brown et al. 1994; Perryman et al. 1997) and 
has an age of ~ 3Myr (Zapatero-Osorio et al. 2002; Oliveira et 
al. 2004). The region has negligible extinction (Bejar et al. 1999; 
Oliveira et al. 2004), and has been extensively studied in the 
optical, X-ray and infrared (e.g., Kenyon et al. 2005; Zapatero- 
Osorio et al. 2002; Jeffries et al. 2006; Franciosini et al. 2006; 
Hernandez et al. 2007; Caballero et al. 2007; Wolk 1996). 
Hernandez et al. (2007) have obtained Spitzer images of a large 
area in <x Ori in the four IRAC bands and with MIPS at 24 pm; 
they find 336 candidate members, of which 66% are class III 
stars and 34% show evidence of disks. 

Accretion rates have been obtained by Gatti et al. (2008) for 
35 objects in <x Ori from the luminosity of the near-IR hydro- 
gen line Pay; they found mass accretion rates lower on aver- 
age than in younger regions. However, their sample was small 
and limited in mass (0.12-0.5 M Q ). In this paper, we present 
the results for a much larger and better characterized sample 
from Hernandez et al. (2007). We measure mass accretion rates 
from the U-band excess emission, which originates in the accre- 
tion shock where accreting matter impacts on the stellar surface 
(Gullbring et al. 1998; Calvet & Gullbring 1998). The corre- 
lation between the U-band excess luminosity and the accretion 
luminosity has been established both empirically (Gullbring et 
al. 1998; Herczeg & Hillenbrand 2008) and theoretically (Calvet 
& Gullbring 1998). The U-band excess is an excellent proxy of 
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the accretion luminosity, which allows obtaining reliable values 
of the mass accretion rate for large samples of stars using little 
observing time when, as in cr Ori, the extinction is negligible. It 
very well complements measurements obtained from other trac- 
ers, such as the IR hydrogen recombination line luminosities (see 
the discussion in Herczeg & Hillenbrand 2008). 

The paper is organized as follows: observations and data re- 
duction are described in Sect. [2] the properties of the sample are 
derived in Sect. [3] in Sect.|4]we discuss the method used to derive 
the accretion properties. The results are discussed in Sect.[5]and 
[6] Three appendices present additional material on: the recom- 
putation of accretion rates in p-Oph with the new distance and 
evolutionary tracks; the accretion properties of BD candidates 
beyond the Spitzer sample; and the properties of transitional and 
evolved disks in our sample. 

2. Observations and data analysis 

2.1. U -band photometry 

We have performed a U-band survey that covers a total field 
of ~1000 arcmin 2 in the cr Orionis cluster. Observations were 
carried out with FORS1 mounted on the UT2 telescope at the 
VLT using uJHIGH filter and were performed in service mode 
during seven nights from October 2008 to March 2009. All the 
nights were photometric, with seeing in the range 0.6-1.5". We 
observed 28 fields (FOV 6'. 8 x 6'. 8) in cr Ori with an expo- 
sure time of 900 seconds each. The distribution of the pointings 
obtained with FORS1 is shown in FigureQ]and the log of the ob- 
servation is provided in TableQ] The long exposure time on each 
field allowed us to reach a U-band limiting magnitude of ~23; 
objects brighter than U~17 mag are saturated. The fields around 
the two brightest stars of the quintuplet system containing cr Ori 
(within 0.2 pc projected distance) had to be excluded because of 
to light contamination. 
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Fig. 1. Fields in the cr Ori cluster that have been imaged with 
FORS1 (solid boxes). Crosses show class III stars (blue in the 
online version of the figure), circles class II, TD and EV ob- 
ject (red in the online version of the figure), as classified by 
Hernandez et al. (2007). 



Table 1. Journal of observations. RA and DEC refer to the point- 
ing center. 



Run dates 


Observed field 


RA 


DEC 


Seeing 


(y-m-d) 




0) 


(°) 


o 


2008-10-05 


field 1 


84.5001 


-2.30238 


0.92 




field 3 


84.7184 


-2.30234 


1.33 


2008-10-06 


field 6 


84.9268 


-2.37146 


1.20 




field 7 


84.9976 


-2.37146 


1.32 




field 8 


84.5002 


-2.47742 


1.51 




field 10 


84.7184 


-2.47742 


1.46 




field 11 


84.8222 


-2.47742 


1.19 


2008-10-08 


field 13 


84.6129 


-2.58167 


0.65 




field 15 


84.8222 


-2.58167 


0.96 




field 16 


84.9268 


-2.58167 


0.95 




filed 18 


84.613 


-2.68734 


1.0 




field 19 


84.7185 


-2.68734 


0.89 




field 20 


84.8222 


-2.68734 


0.99 


2008-12-21 


field 2 


84.613 


-2.30234 


0.81 




field 4 


84.613 


-2.37146 


0.76 




field 5 


84.7184 


-2.37146 


0.72 




field 9 


84.6131 


-2.47742 


0.79 




field 22 


84.613 


-2.79169 


0.74 




field 23 


84.7184 


-2.79169 


1.09 




field 27 


84.7184 


-2.87086 


0.73 




field 28 


84.8221 


-2.89632 


0.68 


2009-01-24 


field 12 


84.9267 


-2.45196 


0.95 




field 21 


84.9267 


-2.68734 


0.96 




field 24 


84.8221 


-2.76621 


1.06 


2009-01-30 


field 17 


84.9976 


-2.55619 


0.99 




field 26 


84.9974 


-2.76621 


0.78 


2009-03-25 


field 14 


84.7653 


-2.57167 


0.70 




field 25 


84.9268 


-2.79169 


0.69 



2.2. Data reduction 

The data were reduced using standard procedures including 
bias subtraction and flat fielding within the IRAF Q package. 
We performed aperture photometry with the PHOT task in the 
APPHOT package, and using noao.digiphot.daophotroutines for 
the photometry extraction. The stellar density was generally low 
enough to make aperture photometry acceptable. The IRAF rou- 
tine MKAPFILE was used to determine and apply aperture cor- 
rections based on ensemble averages of stars in each separate 
frame. Astrometry correction was done to center the telescope 
coordinates. Photometric standard stars from Landolt (1992) and 
Persson et al. (1998) were observed at least once during each 
night and were used to flux-calibrate the images using the task 
PHOTCAL and to set the zero point magnitudes of each ob- 
serving night and for both chips. Aperture photometry was per- 
formed using ten different apertures per image (0.5, 0.6, 0.7, 0.8, 
0.9, 1.0, 1.25, 1.5, 2.0, 3.0 times the average FWHM). The in- 
ner radius of the sky annulus, which allowed us to define the sky 
brightness, was 10 times the average FWHM, while the width of 
the annulus was fixed at 10 pixels. 

The uncertainties on the U-band magnitudes obtained are 
on the order of +0.1 mag. They are dominated by systematic 
errors, the biggest of which is the error on the zero point 
magnitude (about 0.08 mag), which affects all measurements 
in the same manner; a second systematic term, on the order of 
0.01 mag, is due to the color term correction with respect to the 
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u_HIGH filter. Random errors owing to the aperture photometry 
technique used to derive the stellar flux and the sky brightness 
are also very small. 



2.3. U-band variability 

Young pre-main sequence stars are known to be variable, with 
timescale from hours to several days (e.g., Gomez de Castro et 
al. 1998; Hillenbrand et al. 1998; Briceno et al. 2001; Sicilia- 
Aguilar et al. 2005a, 2005b). The U-band variability is probably 
related to variations of the accretion rate, and, although a proper 
study is well outside the scope of this paper, it is interesting to 
estimate how large an effect this is likely to be. 

In our FORS1 data there are 30 objects that lie at the su- 
perposition of two different fields and have therefore been ob- 
served twice, with separations that range between few days and 
few months. Of these, 19 are Class III stars and show no vari- 
ability within the photometric uncertainty. Of the 1 1 Class II ob- 
jects, 5 show no variability, 5 have variations between 0.15 and 
0.4 mag, and one (SO866) has two measurements which differ 
by about 0.8 mag. This is similar to what is observed in other 
star-forming regions (see, e.g., Sicilia-Aguilar et al. 2010 and 
references therein). 

Variability on a much longer time base can be checked by 
comparing our sample to that observed by Wolk (1996), where 
the observations were carried out in January 1996 with the 1 
meter telescope CTIO. The two samples have ~30 class II stars 
in common (see Fig. [2]). The comparison between the U band 
magnitude determinations shows a AU variation of at most 0.5 
magnitudes, not different from what we observe on a shorter 
timescale. 

A difference of 0.5 mag in the measured U-band magnitude 
corresponds to a L acc and M acc difference of a factor of at most 
two (see §4 ). This can be important when discussing individual 
objects. However, if the accretion properties of a large sample of 
stars are considered, it can only cause a moderate spread in the 
accretion values. 
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Fig. 2. U-band magnitude from Wolk (1966) and this work. 
Dashed lines correspond to an interval of +0.5 mag. 



2.4. Spectroscopy 

Among the objects observed with FORS1, six have also optical 
spectroscopy obtained with SARG@TNG. We obtained spectro- 



scopic observations in 2009 at the Telescopio Nazionale Galileo 
(TNG) during three nights from 27 to 29 January. The SARG 
spectrograph attached to the 3.58m telescope was used with 
the 2048x4096 CCD detector (pixel size = 13.5 /mi) and the 
Yellow Grism CD#3 as cross-disperser. This allowed us to cover 
~6200-8000 A wavelength range. We used the slit #1 obtaining 
R -29 000 as spectral resolution. 

The data reduction was performed by using the echelle task 
of the IRAF package following the standard steps of background 
subtraction, division by a flat-field spectrum given by a halogen 
lamp, wavelength calibration using emission lines of a Thorium- 
Argon lamp, and sky subtraction. 

With exposure times of 10-30 minutes we achieved signal- 
to-noise ratios (S /N) in the range 10-15 in the lithium line 
/16708 region, depending on airmass and sky conditions. Each 
star was observed 1-3 times. All the spectra acquired per star 
were shifted in wavelength for the heliocentric correction and 
then co-added obtaining a S /N ratio in the continuum around 
20. 



3. The observed sample 

Table IC21 gives the measured U-band magnitude of all the ob- 
jects in the FORS fields listed as cr Ori members by Hernandez 
et al. (2007) based on optical and near-IR photometry. This sam- 
ple of 187 objects (out of 336 members) is our basic sample and 
will be discussed below. 

The Hernandez et al. (2007) sample spans the mass range 
from ~ 0.06 to 2-3 M , and is practically complete above 0.1 
M . Based on the SED in the IRAC spectral range (from 3.6 fim 
to 8.0 pm ) the Spitzer sources are divided in in class II stars, 
pre-main sequence stars with IR excess typical of optically thick 
disks (classical TTauri stars CTTs or CII), class III stars, with 
typical colors of stellar photospheres (weak-line TTauri stars 
WTTs or CIII), and stars with non-classical disks, in turns di- 
vided in "evolved disks" (EV), with small excess emission at all 
infrared wavelengths, and "transitional disks" (TD), which have 
zero or very low emission in the near infrared but normal excess 
at longer wavelengths. Below we will refer to class II, TD, and 
EV objects as "disk objects". 

Our sample includes 115 class III members, and 72 objects 
with evidence of disks. Among these, there are 54 class II stars, 
4 TD (out of the 7 possibly identified by Hernandez et al. 2007), 
and 14 EV disks. Of the 72 stars with disks, 54 are detected in 
the U-band, 6 are non-detections, and 12 are saturated. Out of 
the 115 class III stars, 83 are detected, 7 are not, and 25 are 
saturated. 

The spatial distribution of the observed sample is shown in 
Fig.[TJ 

3. 1 . Stellar properties 

Spectroscopically determined spectral types exist for a small 
fraction of the cr Ori objects only (Zapatero-Osorio et al. 2002; 
Barrado y Navascues et al. 2003; Muzerolle et al. 2003). 
Therefore, we determine the stellar parameters (effective 
temperature, luminosity, mass, and radius) of all objects in a 
homogeneous way from multi-color photometry. 

Table lC.2l reports for each object broad-band magnitudes col- 
lected from the literature. The optical photometry is from Sherry 
et al. (2004); Zapatero-Osorio et al. (2002); Kenyon et al. (2005); 
Bejar et al. (2001) and Wolk (1996). The JHK magnitudes are 
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taken from the Two Micron All Sky Survey (2MASS) (Cutri et 
al. 2003). The magnitudes in the four channels of the Infrared 
Array Camera (IRAC; 3.8-8.0 pm) and the first channel of the 
Multiband Imaging Photometer for Spitzer (MIPS; 24 pm) are 
from Hernandez et al. (2007). When two or more magnitude 
determinations for the same band were available, we choose, if 
possible, measurements obtained by the same author. 

We derived the effective temperatures of each star by com- 
paring the observed magnitudes to the synthetic colors com- 
puted from the model atmosphere of Baraffe et al. (1998) for 
log g=4.0 and the appropriate filter passband and zero fluxes. 
Luminosities were then computed using the I-band magnitude 
and the bolometric correction for ZAMS stars and the T c yy - 
spectral type correlation of Kenyon & Hartmann (1995) and 
Luhman et al. (2003). We assumed a distance of D=360 pc, and 
negligible extinction in all bands (Brown et al. 1994; Bejar et 
al. 1999). 

We performed two checks on the T f yy estimates. In a first 
test, we compared the values derived from the model atmosphere 
synthetic colors with those obtained by comparing the observed 
colors (V-R), (V-I) and (R-I) to those of ZAMS stars. The dif- 
ferences are within 150 K for 90% of the stars, with a slight 
systematic tendency toward higher values for model atmosphere 
estimates. A second test is given by the comparison of our esti- 
mates of T e ff with the spectroscopic determination derived from 
low- and medium-resolution observations. In a sample of 22 ob- 
jects in the range ~ 3200^4000 K we find differences of + 150 K 
at most. 

The distribution of the stars in the HR diagram is shown in 
Fig. [3] with the evolutionary tracks and isochrones from Baraffe 
et al. (1998). The sample covers the mass range between ~1.2 
to ~0.05 M ; both the lower and the higher mass limit re- 
flect the sensitivity limit of the Spitzer survey (see discussion 
in Hernandez et al. 2007). The median age is about 3 My, with 
a rather large spread, which is similar for Class II and Class III 
objects. This spread also remains when only radial velocity con- 
firmed members are considered (Sacco et al. 2008; Kenyon et 
al. 2005). The issue of the spread in age in cr Ori as well as in 
other young star-forming regions has been discussed in several 
papers (e.g. Hillenbrand 2009 and references therein) and fur- 
ther discussion is beyond the purpose of this paper. Note, how- 
ever, that the error bars may be quite large and affect the age 
estimates significantly. 

The stellar parameters of all Class II objects are summarized 
in Table [2] 



4. Mass accretion rate 

Matter accreting from the disk onto the star, channeled along 
field lines, shocks at the stellar surface. About half of the accre- 
tion luminosity is released with a typical color temperature of 
~ 10 4 K, i.e., much hotter than the stellar photosphere (Hartigan 
et al., 1991, Gullbring & Calvet, 1998). The resulting excess 
emission is clearly detected at short wavelengths, in the U- 
band in particular. It has been shown (Gullbring & Calvet 1998, 
Herczeg & Hillendrand 2008) that the U-band excess luminos- 
ity is an accurate proxy of the accretion luminosity, which can 
be reliably used to measure L acc for T Tauri stars and BDs. 

4.1. Class III sources 

In order to measure the excess of luminosity in U-band for a 
given object, we need to know the measured U-band luminos- 
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Log T eff (K) 

Fig. 3. Location of the observed objects (class III, class II, and 
objects with "non-classical" disks, as labeled) on the HR dia- 
gram. Evolutionary tracks and isochrones are taken from Baraffe 
et al. (1998). Stellar masses and ages are labeled. The horizontal 
bar refers to the error on the effective temperature, the vertical 
bar reflects the error of the bolometric correction of the ZAMS 
stars related to the error on the effective temperatures. 

ity and the expected photospheric U band luminosity for a non- 
accreting star with the same parameters. This is not trivial, as 
young stars tend to have a significant level of chromospheric 
activity that causes continuum emission at short wavelengths; 
this is not related to accretion and must, therefore, be counted as 
"photospheric" contribution for our purpose. 

We take advantage of the large number of Class III stars in 
our sample to define the typical colors of non-accreting young 
stars. Fig. plots the class III (U-J), (U-I), (U-R) and (U-V) 
colors as a function of the effective temperature of the star. For 
each color index, there is a tight correlation with the effective 
temperature; the solid lines show the best fits, the dashed lines 
the 2cr errors, which are ~ 0.5 mag at most. Hereafter, we will 
refer to the +2<x as the photospheric strip. Note that while (U-J) 
increases with T e /f, (U-V) slightly decreases. 

4.2. Disk sources 

To derive the U-band excess emission in class II objects we as- 
sign to each of them the photospheric colors of class III stars 
of the same T e ff> according to the correlations shown in Fig.|4j 
we assume that there is negligible excess emission in the V, R, 
I, and J band and derive the U-band excess from the difference 
between the observed and the photospheric colors. We use (U- 
I), as I-band magnitudes are available for all stars in our sample 
(see Table lC.2"T i. using the relation: 

AUexcess = (U - I) obs - (U - I) pho t, (1) 

where (U - I) bs are the observed (U - I) color, and (U - I) p hot is 
the assigned photospheric color. 

We define as accreting all stars with AU eX cess larger than the 
2cr uncertainties of the photospheric colors, as derived in §4. 1 . 
Class II stars with colors within the class III photospheric strip 
will be considered not-accretors and we can only assign up- 
per limits to the accretion luminosity and mass accretion rate. 
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Fig. 4. Colors of Class III stars vs. effective temperature. Blue 
horizontal lines plot (U-J), red crosses (U-I), green 45deg lines 
(U-R), cyan -45deg lines (U-V), as labeled. The solid lines show 
the best fit, the dashed lines ±2<x. In the text, we will define the 
region between the dashed lines as the photospheric strip. (A 
color version of this figure is available in the online journal.) 



We choose the 2cr uncertainty as good compromise not to loose 
low-accreting stars. After defining the accreting or non-accreting 
stars, we computed the excess flux in the U-band as: 



U,exc 



F ,u x (10 ' 2 - 5 - lO-^-'W-)/"), 



(2) 



where Fqjj is the zero point flux in the U-band. 

Figure [5] shows the comparison between the accretion lumi- 
nosities computed from (U-I) vs those derived from (U-J), (U-R) 
and (U-V), respectively, for all objects with available photome- 
try. The derived L acc are the same within a factor ~ 3 for all ob- 
jects, and for 70% of them the agreement is within a factor 1.5. 
These results support several aspects of our procedure: within 
the above uncertainties: firstly, our assumption that reddening is 
negligible at all wavelengths; secondly, that there is negligible 
excess emission in V, R, I, and J; thirdly, since the photometric 
data are collected from the literature and are not simultaneous, 
that variability is not the major limiting factor in deriving the ac- 
cretion properties (although it may introduce some scatter in the 
measurements). 

We use the (U-I) color to define the accretion rates. Assuming 
an uncertainties in U-band of 0.5 mag we estimate an error on 
the U-band luminosities of a factor 3 at most. 



4.3. Accretion rate 

The U-band luminosity obtained from the excess flux in the 
U-band (Eq. [2]) is converted into total accretion luminosity 
Locc, which is roughly the amount of energy released by gas 
that accretes onto the star, using the approximately linear re- 
lation derived for T Tauri stars (Gullbring & Calvet 1998) 
and brown dwarfs (Calvet & Gullbring 1998; Herczeg & 
Hillendrand 2008): 



iog(^) = io g (^) + i.o. 



(3) 
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Fig. 5. Values of the accretion luminosity derived from (U-J), (U- 
R), (U-V) vs. the value derived from (U-I), as described in the 
text for a sub-sample of 30 objects. The dashed lines correspond 
to ±0.5 in logL acc . 



Knowing M„, R, and the accretion luminosity L fl 
the mass accretion rate M acc with the relation 



R„ 



GM, 



1.25 



L aC c R* 

GM* 



we derive 



(4) 



where G is the universal gravitational constant and the factor 
(1 - f^) -1 ~ 1-25 is estimated by assuming that the accretion 
gas falls onto the star from the truncation radius of the disk (/?,„ 
~5R t ; Gullbring et al. 1998). 



Table [2] reports 72 actual values of M acc . Thirty stars have 
M acc detection. In 30 cases, we can only estimate upper limits to 
Mgcc- Of these, 6 are objects with U-band emission below our de- 
tection limit, while 24 have colors within the photospheric strip 
(§4.2). Eight stars have expected U-band photospheric emission 
brighter than our saturation limit (labeled "sat"), and 4 stars have 
U < 17 mag, but colors with a lower limit than M acc can be safely 
estimated because their photospheric contribution is lower than 
the saturation limit, and the upper limit in the U-band can be 
only due to the accretion process. 

The lowest values of M„ cc we can estimate range from ~ 
10~ u M /y for very low, cold objects to few 10~ 10 M /y for 
solar-mass stars. This trend occurs because the minimum de- 
tectable value of M acc from U-band photometry depends not 
only on the depth of the photometry, but also on the physical 
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Table 2. Accretion properties of the class II objects and objects 
with transitional and evolved disks. 



-3 -2 
Log L acc /L (Pa y) 

Fig. 6. Comparison between L acc computed from Pay (Gatti et 
al. 2008) and L acc obtained from U-band photometry for stars in 
common. The dashed line refers to the same value of L acc , the 
dashed-dotted lines to +0.5 in log L acc . 



properties of the star, and on how well one can estimate the pho- 
tospheric flux. Indeed, this is at present the major limiting factor. 

The uncertainties on individual measurements of L acc and 
M acc are quite large. They come from the combination of photo- 
metric errors and variability, the definition of the class III colors, 
the adopted relation between the U-band excess luminosity and 
the accretion luminosity, and between the latter and M acc and 
the uncertainty on the value of R*/M,(see also Sicilia-Aguilar 
et al. 2010). Some of these uncertainties have been discussed in 
the previous sections, while others, e.g., the differences owing to 
different evolutionary tracks, by other authors (Fang et al. 2009). 
Our estimate is that in general M acc is known with an uncertainty 
of a factor 3-5. 



5. Results 

Detectable values of accretion luminosity and mass accretion 
rate were obtained for 42% of disk objects. Another 42% have 
upper limits in accretion luminosity and mass accretion rate. The 
remaining 16% have lower limits (see §4.3). 

Gatti et al. (2008) compute L acc from the luminosity of the 
hydrogen recombination line Pay, following the procedure de- 
scribed by Natta et al. (2006). 

Figure [6] compares the values of L acc derived with the two 
different methods (U-band excess emission and Pay recombi- 
nation line respectively) for the 12 stars we have in common; it 
agrees well. Another nine stars of the Gatti et al. (2008) sample 
are class II objects and were observed with Spitzer. We could 
not observe these stars either because they are outside the FORS 
fields, or because they are close to the brighter stars of the cr Ori 
quintuplet system. We added these nine stars to our sample. In 
Table. [3] the properties of the sub-sample derived from Gatti 
et al. (2008) are listed. For these stars we recomputed the 
stellar properties adopting the Baraffe et al. (1998) evolutionary 
tracks., in the same way as for the present sample. 

Figures [7J - [10] show the relations between the accretion 
properties and different physical and morphological properties 
of the observed sample. 

In Fig. [7J we specifically plot the accretion luminosity as 
a function of the stellar luminosity for all class II objects; 
transitional and evolved disks are shown by different symbols. 
The figure shows that very few stars have L acc larger than 0. 1 
L star , and that most of them have values well below this limit. 



Object 


class 


M» 


L, 


Tcff 


logL acc 


logM acc 






(M ) 


(L ) 


(K) 


(L Q ) 


(M yr-') 


S0299 


TD 


0.35 


0.14 


3400 


-2.93 


-9.82 


S0341 


II 


0.60 


0.39 


3600 


< -2.26 


< -9.32 


S0397 


II 


0.30 


0.18 


3300 


-2.73 


-9.48 


S0435 


II 


0.20 


0.09 


3200 


< -3.38 


< -10.07 


S0444 


EV 


0.35 


0.16 


3400 


< -2.90 


< -9.78 


S0462 


II 


0.20 


0.15 


3200 


-3.08 


-9.67 


S0482 


II 


0.20 


0.05 


3200 


-3.08 


-9.91 


S0485 


II 


0.25 


0.04 


3300 


-2.39 


-9.37 


SO490 


II 


0.15 


0.07 


3100 


-2.73 


-9.32 


SO500 


11 


0.10 


0.02 


3000 


-3.27 


-9.92 


S0514 


II 


0.13 


0.04 


3100 


< -3.88 


< -10.56 


S0518 


II 


1.00 


0.41 


4000 


sat 


sat 


SO520 


II 


0.30 


0.15 


3300 


< -3.02 


< -9.82 


S0537 


11 


0.07 


0.02 


3000 


< -4.32 


< -10.83 


SO540 


11 


1.00 


0.43 


4000 


sat 


sat 


S0562 


II 


0.40 


0.22 


3400 


-1.80 


-8.66 


S0563 


II 


0.40 


0.51 


3300 


-2.30 


-8.95 


S0583 


II 


0.75 


1.52 


3600 


sat 


sat 


S0587 


EV 


0.30 


0.27 


3300 


< -2.78 


< -9.45 


S0598 


II 


0.20 


0.11 


3200 


-3.09 


-9.75 


S0615 


EV 


0.60 


0.69 


3500 


sat 


sat 


S0638 


EV 


1.00 


0.46 


3900 


sat 


sat 


S0646 


11 


0.30 


0.12 


3300 


-2.05 


-8.89 


S0657 


II 


0.06 


0.02 


2900 


< -4.49 


< -10.91 


S0662 


II 


0.90 


0.41 


3900 


< -2.06 


< -9.26 


S0663 


11 


0.20 


0.14 


3200 


< -3.20 


< -9.82 


S0674 


II 


0.30 


0.14 


3300 


< -3.06 


< -9.88 


S0697 


II 


0.60 


0.63 


3500 


-1.92 


-8.76 


SO700 


EV 


0.08 


0.01 


3000 


< -4.46 


< -11.11 


SO710 


II 


0.50 


0.30 


3500 


-2.28 


-9.19 


S0728 


EV 


0.20 


0.25 


3100 


< -3.08 


< -9.53 


S0736 


II 


0.80 


1.73 


3700 


sat 


sat 


S0738 


II 


0.11 


0.03 


3100 


< -4.04 


< -10.7 


S0739 


II 


0.10 


0.04 


3000 


< -4.03 


< -10.54 


SO750 


II 


0.09 


0.03 


3000 


-3.51 


-10.03 


S0759 


EV 


0.30 


0.19 


3300 


< -2.97 


< -9.67 


S0762 


II 


0.13 


0.05 


3100 


< -3.82 


< -10.44 


S0774 


II 


1.00 


0.47 


4000 


< -1.93 


< -9.17 


S0818 


TD 


0.70 


0.21 


3700 


-2.26 


-9.45 


S0827 


II 


0.35 


0.06 


3400 


-2.74 


-9.82 


S0844 


II 


0.60 


0.42 


3600 


-2.01 


-8.95 


S0848 


II 


0.20 


0.03 


3300 


-2.85 


-9.86 


S0859 


II 


0.35 


0.17 


3400 


-2.75 


-9.61 


S0865 


II 


0.35 


0.13 


3400 


-2.59 


-9.50 


SO866 


II 


0.20 


0.05 


3200 


-3.28 


-10.12 


S0897 


TD 


0.80 


0.50 


3700 


> -2.07 


> -9.13 


SO905 


EV 


0.60 


0.35 


3600 


< -2.39 


< -9.38 


SO908 


EV 


0.20 


0.10 


3200 


-2.68 


-9.37 


S0917 


EV 


0.20 


0.03 


3300 


< -3.73 


< -10.71 


S0927 


II 


0.70 


0.29 


3700 


< -2.41 


< -9.53 


S0936 


II 


0.08 


0.01 


2900 


< -4.58 


< -11.19 


S0967 


II 


0.25 


0.08 


3300 


< -3.39 


< -10.25 


S0981 


EV 


1.3 


1.17 


4200 


sat 


sat 


S0984 


II 


0.90 


0.53 


3800 


< -2.12 


< -9.24 


SO1009 


EV 


0.15 


0.02 


3200 


-1.97 


-8.89 


SO1036 


II 


0.75 


0.45 


3700 


< -2.17 


< -9.22 


SO1050 


II 


0.20 


0.08 


3200 


< -3.49 


< -10.22 


SO1057 


EV 


0.15 


0.08 


3100 


< -3.67 


< -10.25 


SO1059 


II 


0.06 


0.01 


2900 


< -4.33 


< -10.86 


SO1075 


II 


0.17 


0.13 


3100 


> -1.66 


> -8.18 


SOH56 


II 


1.00 


0.63 


3900 


sat 


sat 


S01182 


II 


0.20 


0.10 


3200 


-3.12 


-9.79 


S01193 


II 


0.11 


0.03 


3100 


-3.32 


-10.01 


SO1230 


II 


0.20 


0.08 


3200 


-3.32 


-10.06 


SO1260 


II 


0.35 


0.13 


3400 


-1.68 


-8.60 


S01266 


II 


0.16 


0.07 


3100 


-3.47 


-10.09 


S01268 


TD 


0.09 


0.02 


2900 


< -4.49 


< -11.06 


S01274 


11 


0.70 


0.64 


3600 


> -2.15 


> -9.07 


S01323 


EV 


0.20 


0.07 


3200 


< -3.51 


< -10.25 


S01327 


II 


0.30 


0.20 


3300 


-2.66 


-9.39 


S01361 


II 


0.55 


0.53 


3500 


> -1.82 


> -8.65 


S01362 


II 


0.15 


0.09 


3100 


-3.07 


-9.61 
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Table 3. Accretion properties of the class II objects taken from 
Gatti et al. (2008). The masses and luminosities of the listed stars 
have been recomputed according the Baraffe et al. (1998) evolu- 
tionary tracks. 



Object 


class 


M, 


L, 


T e ff 


logL acc 


logM acc 






(M ) 


(Lq) 


(K) 


(L ) 


(Moyr 1 ) 


S0451 


II 


0.30 


0.12 


3350 


-2.20 


-9.10 


S0682 


11 


0.60 


0.30 


3575 


< -2.50 


< -9.30 


S0694 


II 


0.15 


0.10 


3075 


< -3.20 


< -9.80 


S0723 


II 


0.25 


0.10 


3250 


-2.00 


-8.90 


S0726 


11 


0.60 


0.42 


3575 


-1.80 


-8.60 


S0733 


II 


0.45 


0.36 


3425 


-2.00 


-8.70 


S0871 


II 


0.45 


0.10 


3500 


< -3.20 


< -10.30 


SOH52 


11 


0.60 


0.30 


3575 


< -2.40 


< -9.30 


SO 1248 


n 


0.13 


0.13 


3000 


< -2.80 


< -9.30 




Fig. 7. Accretion luminosity as a function of the stellar luminos- 
ity for class II, EV, and TD disks, as labeled. The dashed lines 
show L acc IL star = 0.001, 0.01, and 0.1, respectively. The crosses 
surronded by a circle are the class II stars included here from the 
Gatti et al. (2008) sample, and listed in Table [3] Different col- 
ors for lower and upper limits can be distinguished in the online 
version of the journal. 



The accretion luminosities for the detected sources range mainly 
from ~ 10~ 2 L Q to ~ 10~ 4 L o . For any given L„ there is a large 
range of measured L acc that does not seem to vary significantly 
with L„ . 

Figure [8] shows the mass accretion rate of class II stars as 
function of M,. The data show a clear trend of increasing M acc 
with increasing M». Including upper and lower limits as actual 
detections, we find M acc oc M\- b±aA with ASURV (Astronomy 
Survival Analysis Package, Feigelson & Nelson 1985). The 
trend is confirmed using two different methods (the EM algo- 
rithm and the BJ algorithm) within the ASURV package. The 
slope became flatter when we excluded the upper and lower lim- 
its in the analysis, but remained still within the uncertainties. 
From this plot we can clearly also see the large spread in M acc 
(about two orders of magnitude) for any value of M*. 

In Fig [9] we plot the mass accretion rates versus the age of 
the stars. 

Finally, in Fig [10] we plot the mass accretion rates as func- 
tion of the projected distance from the central and bright 09.5 
star cr Ori. We do not find any correlation between the mass 
accretion rates or the disk morphologies and the distance from 
cr Ori; it seems that the vicinity of the 09.5 star does not affect 
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Fig. 8. Mass accretion rate as function of stellar mass. Symbols 
as in Fig. [7] (A color version of this figure is available in the 
online journal.) 




6.2 6.4 6.6 6.8 
Log (Age) (My) 

Fig. 9. Mass accretion rate as function of the age of the stars. 
Symbols as in Fig. [7] 



the disk properties significantly (see also Hernandez et al. 2007). 
Rigliaco et al. (2009) studied in detail one object (S0587), which 
has a projected distance of 0.3 pc; based on literature U-band 
photometry and optical spectroscopy, they proposed that its disk 
is in the process of being photoevaporated, either by cr Ori or 
by its own central star. There are six additional objects with a 
projected distance of<0.3 pc from cr Ori with a M acc range of 
more than a factor of ten, and it would be interesting to obtain 
high-resolution optical spectra to study wind diagnostics such as 
the optical forbidden lines of [SII], [Nil] and [OI]. 

We caveat that in both Fig. [9] and Fig. [10] we did not divide 
stars in mass bins (which would result in too low number statis- 
tics); thus the relationship of M acc vs. M* seen in Fig. [8] might 
mask possible trends between M acc and age or M acc and the pro- 
jected distance from cr Ori. 



6. Discussion 

6.1. M occ as function of Af« and time 

The distribution of mass accretion rates in a star-forming region 
traces the physical processes that control disk formation and evo- 
lution over the lifetime of the region, as well as the initial con- 
ditions, i.e., the mass and angular momentum distribution of the 
molecular cores from which the stars form. It is a snapshot in 
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Fig. 10. Mass accretion rates versus projected distance from the 
bright stars cr Ori in parsec. (A color version of this figure is 
available in the online journal.) 



time, which needs to be compared to models that follow disk for- 
mation and evolution up to the age of the region (e.g., Dullemond 
et al. 2006; Vorobyov & Basu 2008, 2009). 

It has been known for some years that mass accretion rates 
increase on average with the mass of the central object. In <x Ori, 
a logarithmic linear correlation over the whole mass range, 
M acc cc Ml' 6±0A (Sec. 01, is similar within the uncertainties to 
what has found in Taurus (Calvet et al. 2004) and Ophiuchus 
(Natta et al. 2006), but flatter than in L1630N and L1641 (Fang 
et al. 2009). However, in cr Ori the trend seems to be flatter for 
higher mass stars than for lower mass objects, as also noted by 
Vorobyov & Basu (2008) in a compilation of all known M acc 
values in various star-forming regions. This is shown in Fig.fTTl 
which plots median values of M acc as function of M, for Class 
II and TD disks. All M* intervals above 0.1 M Q contain roughly 
the same number of stars (12-20). The values plotted in Fig.Qj] 
treat upper and lower limits as actual detections; if we make the 
alternative assumption that lower (upper) limits are all smaller 
(larger) than the lowest (highest) measured value in the bin, the 
changes are very small; however, the large number of limits and 
the small number of objects in each bin make it meaningless 
to derive values for the upper and lower quartiles, to charac- 
terize the spread of M acc seen in Fig|8] The mass interval be- 
low 0. 1 M contains only two detections and six upper limits to 
Macc, and the value in Fig.[TT](Log M flCC =-10.8 M /y) should be 
considered as an upper limit to the median. The deficit of rel- 
atively strong accretors among the <x Ori BDs is confirmed by 
the analysis of a much larger (40 objects, see Appendix B), op- 
tically selected sample of very low-mass stars and BDs (Lodieu 
et al. 2009), included in our U-band survey, which contains only 
five objects with Log M acc > -10.8 M /y, two of which are also 
in the Hernandez et al. (2007) Spitzer sample. 

The flattening of the M acc vs. M, relation at higher masses is 
very clear. Vorobyov & Basu (2009) compute numerical models 
of the collapse of a distribution of prestellar cores that include 
the formation and evolution of circumstellar disks and follow 
it for 3 My, roughly the age of <x Ori. Their models predict a 
flattening of the M acc vs. M, relation for stellar masses higher 
than about 0.3 M , due to the effect of gravitationally induced 
torques in the early stages of the evolution after the formation 
of the central star. These gravitational instabilities have little ef- 
fect on lower-mass objects, where viscous evolution dominates 
at all times. Although the M acc values predicted by Vorobyov 



& Basu (2009) are somewhat higher than the observations, the 
cr Ori results definitely support their models and the importance 
of self-gravity in the early evolution of more massive disks, al- 
ready suggested by, e.g., Hartmann et al. (2006). 

The Vorobyov & Basu (2009) models include only viscous 
evolution and gravitational instability; other physical processes 
may occur during cr Ori lifetime, such as photoevaporation and 
planet formation, leading to disk dissipation on shorter time 
scales. This may be a selective process, if, as indicated by 
the statistics of IR-excess emission for stars of different mass 
in different star-forming regions, disk dissipation occurs faster 
in more massive stars. In cr Ori, for example, Hernandez et 
al. (2007) estimate a fraction of objects with disks that increases 
from ~ 10% for Herbig Ae/Be stars to ~ 35% for T Tauri stars 
and BD candidates. If so, the comparison of the observations 
with the model predictions needs to be taken with care. 

A way to investigate the relative importance of these dif- 
ferent processes is to compare the statistical properties of the 
distribution of M acc on M* for regions of different age. This 
approach is limited at the moment because to the best of our 
knowledge, there are only two other suitable samples: p-Oph 
(Natta et al. 2006) and Tr 37 (Sicilia-Aguilar et al. 2010). In 
other cases, no mass accretion rates are available for complete 
sample of Class II objects (as in Taurus), or upper limits to M acc 
are not provided (as in the two Orion regions studied by Fang et 
al. 2009), which makes statistical studies very difficult. 

We have computed median values of M acc in p-Oph from 
the results of Natta et al. (2006), which we revised to take into 
account the new estimates of the distance that were recently pub- 
lished (see Appendix A for details). The p-Oph sample covers a 
mass range between ~ 0.03 and 3 M ; the median values of M acc 
are shown in Fig. Q~T] the figure also shows the results for the 
older (about 4 My), more distant region Trumpler 37, for which 
Sicilia-Aguilar et al. (2010) provide values of M acc derived from 
U-band photometry for a Sp/fzer-selected sample of stars. The 
stars cover the mass range 0.4-1.6 M , based on the Siess et 
al. (2000) evolutionary tracks. 

For M, roughly larger than 0.2 M , the three regions have 
similar values of the median M acc within the uncertainties in 
spite of their difference in age (more than a factor of 3). If 
the disk properties at a very early stage were the same, this 
would imply a slower time evolution of M acc than predicted by 
disk models (see, e.g., Hartmann et al. 1998; Dullemond et al. 
2006; Vorobyov & Basu 2009) and confirms the result of Sicilia- 
Aguilar et al. (2010), based on a small sample of stars for which 
individual ages could be estimated. The difference between M acc 
medians increases as M* decreases, and becomes very large (fac- 
tor of 10 at least) for M* < 0.1 M, between p-Oph and cr Ori. 
This is consistent with the predictions of viscous models if the 
original disk properties in the two regions are the same and no 
disk around these very low-mass stars dissipates within, e.g., 3 
My. 

If indeed the apparent slow time evolution of solar mass stars 
is a result of the continuous loss of the less massive, for the lower 
accreting disks the two different effects (namely the decrease of 
M acc with time, and dissipation of disks once they fall below 
a critical M acc value) roughly compensate, so that median M acc 
changes little. If that is true, one should find that the lowest val- 
ues of M acc do not vary with time, because they are fixed by disk 
dissipation, while the highest values decrease as expected from 
viscous evolution. There is a hint that this is indeed the case (seel 
also L1641 and L1630; Fang et al. 2009), but the statistics is poor 
and the definition of the upper envelope to the M acc values too 
uncertain for the moment. 
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Fig. 11. Median values of M acc as function of M,. Red lines 
(with a central square) refer to cr Ori Class II and TD disks; the 
lowest mass bin should be considered as an upper limit only (see 
text). Blue lines (crosses) show the distribution for p-Oph (data 
from Natta et al. 2006, see Appendix A); green lines (circles) for 
Tr 37 (data from Sicilia-Aguilar et al. 2010). (A color version of 
this figure is available in the online journal.) 



6.2. Mass accretion rates vs disks morphologies 

Fig. Q~2] plots the mass accretion rates for the cr Ori sample 
as function of the excess over the photospheric emission for 
the four IRAC bands and for the 24 /mi MIPS band. The ex- 
cess emission is defined as the difference between the observed 
(Spitzer-T) color and that of Class III stars of the same effective 
temperature, assuming that no excess is present in the I band. 
In each panel, the horizontal bars show the median of the ex- 
cess distribution for Class II disks and the lower and upper first 
quartile; for comparison, we also show the lower and upper first 
quartile of the excess for the Taurus median SED of classical 
disks (D'Alessio et al. 2006). 

The median of the excess emission distributions is lower in 
cr Ori than in Taurus, as expected if on average cr Ori stars are 
older than Taurus objects, suggesting a higher degree of grain 
settling in older regions (e.g., Hernandez et al. 2007). Note that 
the cr Ori sample includes lower-mass objects than the Taurus 
one, and that disk models predict lower excess fluxes for BDs 
in this range of wavelengths; however, we do not find any cor- 
relation of the measured excess with the mass of the star, and 
we tend to exclude that the difference between the two regions 
is due to the different mass range only. The figure also shows 
the location of the 1 1 EV disks with a determination of M acc and 
the 4 TD; for six of them we also show complete SEDs and Ho- 
profiles in Appendix C. 

There is no statistically significant correlation of M acc with 
the excess emission. However, there are a few aspects of these 
plots that can be understood if, as expected in viscous mod- 
els, M acc traces the surface density of the inner disk. Objects 
with M acc ~ 10" 10 - 10~ 9 M /y are distributed over the whole 
range of excess values. This is expected in optically thick disks, 
where the observed emission depends on inclination, inner ra- 
dius and degree of flaring, but not on the actual disk surface 
density. Objects with M acc <10~ 10 M Q /y tend to have very low 
excess in all Spitzer bands. We think that most of them have 
in fact optically thin disks, which are characterized by much 
lower emission, roughly proportional to the surface density, i.e., 
10" M G /y is a reasonable threshold 



for the transition from optically thick to optically thin (inner) 
disks (D'Alessio et al. 2006). 

A last point to note is that disks with M flC( >10~ 9 M Q /y have 
all large excess emission (in the upper quartile of the distribu- 
tion), and indeed the upper envelope of the M acc distribution 
seems to correlate with the amount of excess emission. Given 
the small number of objects in the high range of M acc , the sig- 
nificance of this is unclear. If true, it would suggest a very inter- 
esting relation between M acc and grain growth and settling, i.e., 
processes that can change the grain opacity and the disk flaring. 

The potential of plots like Fig.Q~2]in constraining disk model 
parameters should be exploited further. However, this is well be- 
yond the scope of this paper. 

Fig. [12] also shows the location of the 11 EV disks with an 
M acc determination. Nine of them have upper limits to M acc , 
some well below 10 M G /y, and very small excess emission, 
generally below the lower quartile of the Class II distribution. 
From the upper limits to M acc , it is likely that at least half of them 
are indeed optically thin disks. There are two exceptions, one 
(SO908) is probably a missclassified Class II object, as it shows 
significant excess emission in all bands and broad Ha emission, 
consistent with its measured accretion rate (M acc - 4 x 10~ 10 
M Q /y) (see Appendix C). The other EV object (SO1009) has one 
of the highest accretion rates in our sample (1.2x1 0~ 9 M Q /y) and 
no detectable excess emission in all bands; indeed, its classifica- 
tion as a disk object is dubious. No additional data are available 
in the literature, and we may have detected a strong chromo- 
spheric flare. It would be interesting to monitor this object fur- 
ther. 

The location of the four TD on Fig. [12] is also shown. They 
are scattered through the plot: two (S0818 and S0897) are con- 
sistent (both in M acc and excess emission, see also Appendix C) 
with optically thick inner disks. One (S01268) has very likely 
a very optically thin inner disk, with a very low upper limit to 
M aC c{~ 9 x 10~ 12 M /y) and small excess emission even at 24 
pm. The fourth TD object (S0299) has an accretion rate typical 
for its mass, 24 pm excess as Class II of similar accretion rate, 
but negligible excess in the IRAC bands. The large spread of 
properties of TD confirms the analysis of Muzerolle et al. (2010) 
and Sicilia-Aguilar et al. (2010) and their conclusions that very 
likely the "transitional" SEDs trace a variety of different physi- 
cal situations. 



7. Conclusions 

We reported the results of a U-band survey with FORS1/VLT 
of a large area in the cr Ori star-forming region. We combined 
the U-band photometry with literature results to compute accre- 
tion luminosity and mass accretion rates from the U-band excess 
emission for all objects detected by Spitzer in the FORS1 field 
and classified by Hernandez et al. (2007) as likely members of 
the cluster. In total, there are 72 objects with evidence of a disk 
from near- and mid-IR photometry and 115 class III (diskless) 
stars. Among the disk objects, four (out of the seven identified by 
Hernandez et al. 2007) are transitional disks, and 14 are evolved 
disks. We derived the photospheric parameters of all stars from 
the existing V, R, I, and J photometry and used the U-A colors of 
class III as templates for the photospheric and possible chromo- 
spheric emission. Our final sample, for which we provide esti- 
mates of the mass accretion rates, contains 58 Class II (49 class 
II stars for which we derive the accretion properties from the U- 
band excess emission, and nine stars with accretion properties 
from literature, for which we checked the consistency with our 
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Fig. 12. M ace vs excess emission (in mag) in the four IRAC 
Spitzer bands and for the 24pm MIPS band. The median excess 
emission and first quartiles are shown in each panel by the red 
(thick) horizontal lines; for comparison, we also plot the median 
Taurus class II SED (black (thin) line; D'Alessio et al. 2006). 
Symbols as in Fig. [7] (A color version of this figure is available 
in the online journal.) 

results), four TD and 1 1 EV disks, over a mass range of between 
~ 0.06 and ~ 1.2M . 



We analyzed the behavior of M acc as function of mass and 
age of the individual stars, of the properties of the IR SED, and 
of the distance from the bright star cr Ori. 

There is no correlation of M acc with the distance from cr Ori, 
confirming that the effect of the 09.5 star on its surroundings is 
not strong (see also Hernandez et al. 2007); however, our sample 
does not include objects with projected distances smaller than 
0.1 pc. 

We find a strong relation between M acc and M», with a very 
large spread of M acc values for any given M„ similar to other 
star-forming regions (Calvet et al. 2004; Natta et al. 2006; Fang 
et al. 2009; Sicilia-Aguilar et al. 2010). If fitted with a linear 
correlation, the slope is 1.6 + 0.4. As noted by Vorobyov & 
Basu (2009) for a compilation of accretion rates in different star- 
forming regions, a linear fit is not the best description of the 
data. We computed median M acc values and showed that the 
relation between M acc and M» is flatter at higher masses and 
steepens significantly for very low-mass stars and BDs. Such a 
trend is predicted by models of core collapse and disk evolution 
that include viscosity and gravitational instabilities (Vorobyov & 
Basu 2009), which control the evolution of more massive disks. 

These models follow the disk evolution to the age of 
cr Ori, under the assumption that disks evolve only by accret- 
ing onto the star. However, other processes such as photoevapo- 
ration and/or planet formation, may cause disk dissipation (e.g. 
Hollenbach 2000; Dullemond et al. 2007). This may be a se- 
lective effect, which affects higher mass stars more than lower 
mass stars, changing the observed dependence of the M acc dis- 
tribution with M t ■ We have compared the M acc -M* distribution 
in cr Ori to that of the two other star-forming regions in the litera- 
ture for which an IR-selected sample is available, namely p-Oph 
(Natta et al. 2006 and Appendix A) and Tr 37 (Sicilia-Aguilar et 
al. 2010). The comparison indicates that the median M acc values 
for higher M, are closer than predicted by simple viscous mod- 
els, suggesting that selective disk dissipation may be important. 
However, we note that the significance of this comparison is not 
very strong, because the results for p-Oph are uncertain and the 
number of objects in each mass bin is not large. Moreover, the 
number and distribution of upper limits to M acc affects the deter- 
mination of median values in some mass ranges and prevents us 
from deriving upper and lower quartiles of the distributions. 

The behavior of M acc as function of the excess emission in 
the Spitzer bands suggests that at M acc ~ 10~ 10 M /y the (inner) 
disks change from optically thin to optically thick. Objects with 
Macc m me range ~ IO" 10 - 10~ 9 M Q /y span the whole range 
of observed excesses, from very low to very large. Objects with 
M acc ^10 -9 M /y (the largest values in the sample) all have large 
IR excess. Viscous disk models (e.g., D'Alessio et al. 2006) pre- 
dict that for M flcc =10~ 10 the surface density at 1 AU from a T 
Tauri star will be on the order of 1-10 g cm 2 , depending on the 
grain properties and dust settling. The emission of these disks 
in the Spitzer bands will be optically thin, unless grains are sub- 
micron size. The trend we tentatively observe among optically 
thin and optically thin disks, which needs to be confirmed in 
larger samples, may indicate a link between the mass accretion 
rate and the grain properties, which in turn control the disk ge- 
ometry,a connection that is worth to be further explored. 

The four TD stars included in our sample seem to cover a 
variety of properties, and only one of them has a mass accretion 
rate as Class II of similar mass, and negligible excess emission to 
wavelengths > 8//m. The other three could not be distinguished 
from Class II objects in the M acc vs. excess emission plots. We 
can only agree with the conclusions that the "transitional" prop- 
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Fig.A.l. Mass accretion rate versus stellar masses for the star 
of the p-Oph sample. The stellar parameters and the accre- 
tion properties were recalculated assuming a distance of 130 pc 
(Lombardi et al.,2008) instead of 160 pc assumed by Natta et 
al. (2006), and 1 My age. 



erties of the SEDs are likely caused by a variety of different prop- 
erties (Muzerolle et al. 2010; Sicilia-Aguilar et al. 2010). 
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Appendix A: p Ophiucus 

The p-Oph sample is particularly interesting, because it is sim- 
ilar to the cr Ori one in being an IR-selected sample of class 
II, complete to a limiting mass of about 0.05 M Q (Bontemps 
et al. 2001). Natta et al. (2006) computed mass accretion rates 
from the luminosity of Pa/? assuming a distance of 160 pc; since 
no extensive spectral type determinations were available, they 
followed the method outlined by Bontemps et al. (2001) and de- 
rived the stellar parameters assuming coeval star formation at 0.5 
My and the evolutionary tracks of D'Antona & Mazzitelli (1998; 
the isochrone method). They found that the mass accretion rates 



,1.8±0.2 



over a mass tn- 



could be fitted by a linear relation oc M- 
terval 0.03-3 M Q , with a very large spread for any given M*. 
However, the Natta et al. (2006) results need to be reconsidered, 
since new measurements of the p-Oph distance yield consider- 
ably low values, 120-130 pc (Lombardi et al. 2008; Loinard et 
al. 2008; Snow et al. 2008). The corresponding decrease in lu- 
minosity implies an older age for the region. We have redeter- 
mined stellar parameters and mass accretion rates for all objects 
in Natta et al. (2006), adopting a distance of 130 pc and evolu- 
tionary tracks of Baraffe et al. (1998) for 1 My. The new values 
of M acc are somewhat lower than previous ones, and M» higher, 
especially for higher masses (see Fig. lA.U . 

As a consequence, if fitted with a single power-law, the cor- 
relation between M acc and M* is flatter, with a slope of 1 .3 + 0.2. 
The new distance, the older age and the different evolutionary 
tracks contribute to this result. In particular, the dependence of 
the M acc - M* relation on the adopted evolutionary tracks is well 
known (see Fang et al. 2009) and is particularly strong in p-Oph, 
given the method used to determine the stellar parameters. 



Appendix B: BD and very low-mass stars from the 
Lodieu et al. photometric survey 

An independent sample of very low-mass stars and brown dwarfs 
has been selected from the list of cr Ori members and candidate 
members of Lodieu et al. (2009). We applied a first selection cri- 
terium on the z,(z-J) diagram computed T e ff by comparing ob- 
served z, Y, J colors to synthetic ones from the theoretical models 
of Baraffe et al. (1998) for logg=4.0 and luminosities from the 
observed J mag and model-predicted bolometric corrections. We 
then performed a further selection based on the location of the 
objects on the HR diagram, excluding all stars with M»>0.13 
M Q . Our sample of candidate young very low-mass stars and 
BDs in cr Ori is then formed by 80 objects, 40 of which were 
included in the U-band FORS 1 survey. 

Spitzer IRAC detections exist for 21/40 objects; three 
of them are uncertain members according to Hernandez et 
al. (2007). Of the 21, six are classified as class II, one is a tran- 
sitional disk, one an evolved disk and 13 are class III objects. 
The 18 confirmed members are included in the sample analyzed 
in the main text of the paper. Note that the determination of the 
stellar parameters, T e ff in particular is performed using different 
photometric bands with respect to the bands used in the main 
text; the difference in T c yy are in general within the uncertain- 
ties discussed in §3.1; however, some of the objects in Table 2, 
although included in the Lodieu sample, were not selected with 
the criteria applied here. 

Of the 40 objects in our sample, 20 have been detected in 
the U-band, while for the other 20 we have upper limits only. 
We derived the accretion luminosity and mass accretion rates 
as in §4.3. The calibration of the photospheric colors (U-J) and 
(U-I) as function of T e ff using our U-band photometry of Class 
III objects extends to T?//~ 2900 K (Sec.4.1); we compare it to 
synthetic colors from the Baraffe et al. (1998) model atmosphere 
to extend the relations to lower T e ff- The Class III colors agree 
well with the models; the results are shown in Fig. IB. II 

Fig. IB. 21 shows M acc vs. M*; we find that only five objects 
(all with M„ 0.06 M ) have M acc > -10.8, the median in §6.1. of 
these, two are class II sources (SO500 and S0848, also in Table 
2), one is a class III (S0641, possibly misclassified), three have 
no Spitzer detections. No other object with higher values of M acc 
is detected. 



Appendix C: SEDs and H Q of EV and TD objects 

In this section, we show SEDs (Fig. IC.U and Ha profiles 
(Fig. IC.2I > of a subset of two TD and four EV stars observed 
with SARG@TNG (Sect. 2.4) and Giraffe (Sacco et al. 2008). 
A summary of the Ha properties is given in Table IC.U where 
10% Ho- represents the width of Ha at 10% of the line's peak in- 
tensity, and pEW is the pseudo-equivalent width. Below we will 
briefly comment on each object. 

S0587. This EV disk was extensively studied by Rigliaco et 
al. (2009). It exhibits modest excess emission (well below the 
lower quartiles of the distributions) in the IRAC bands and at 24 
yum. It shows a symmetric and narrow Ha emission profile with 
the peak close to the line center (type I profile following the clas- 
sification of Reipurth et al. (1996)). Based on the available U- 
band photometry (Wolk 1996), the narrow Ha and the strength 
and profiles of the [S n] and [N n] forbidden lines, Rigliaco et al. 
(2009) suggested that the disk was being photoevaporated and 
that the forbidden lines were coming from the photoevaporation 
flow, possibly driven and certainly illuminated by the star cr Ori. 
A crucial ingredient of this model was the high ratio between the 
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Fig. B.l. U-A colors as function of T e /f. The red circles are 
objects with measured U-band fluxes, arrows are objects with 
U-band upper limits. The squares are model-predicted colors 
for different gravity, from 5.5 (black squares, top) to 3.5 (blue 
squares, lowest). The solid lines show the best-fit relations for 
Class III derived in Sec.4.1, extrapolated to lower T e ff> dashed 
lines are +2cr. We will consider as accretors objects with col- 
ors below the photospheric strip: five BDs have clear evidence 
of U-band excess (note that for two we only have two colors), 
one only marginal evidence (from (U-I) and (U-J), while (U-Z) 
is photospheric). Hereafter we mark the upper limit as dashes for 
clarity. (A color version of this figure is available in the online 
journal.) 
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Log M»/M G 

Fig. B.2. Mass accretion rate as function of M*. Dots are class 
II objects with measured M acc , the cross is a class III stars with 
measured M acc , and diamonds refer to object with no Spitzer 
data. Arrows are objects with upper limits (both U-band detec- 
tions and non-detections). Colors (only in the on-line version) 
indicate the SED Class : red for class II, blue for class III; black 
for objects with no Spitzer data. (A color version of this figure is 
available in the online journal.) 



mass-loss and the mass-accretion rate. Our results confirm this 
interpretation. Although relatively bright in U, the object does 
not have a measurable U-band excess (i.e., it lies inside the pho- 
tospheric strip of Fig. |4]i and we estimate an upper limit to log 
Macc of -9.45 M Q /y, consistent with the value -9.52 M Q /y de- 
rived in Rigliaco et al. (2009). 

S0615. This is a relatively massive, luminous star. Its SED 
is typical of a flat disk up to 24 pm. We cannot measure M acc 
from the U-band, which is already saturated by the photospheric 
emission alone, but we can estimate a value of log M acc =-8.25 
M /y from the U=15.96 mag measurement of Wolk (1996). The 
Ha has a complex profile, with broad wings, deep redshifted and 
blueshifted absorption and a narrow, slightly redshifted emission 
in the center. Similar profiles are observed for higher numbers 
of the Balmer series ("YY Orionis like profiles"; Walker, 1972) 
and are associated with extensive infall and outflow rates, consis- 
tent with the very high value of M acc . Unfortunately, the SARG 
S0615 spectrum is rather noisy. 

S0759. Classified as EV disk star by Hernandez et 
al. (2007). The SED has negligible excess up to 8 pm, but a sig- 
nificant one at 24 ^m, not very different from the SED of some 
TD objects. The Ha is rather narrow and symmetric, with a mod- 
erate red/blue asymmetry (type I profile, Reipurt et al. (1996)). 
The 10% Ha width of 157 + 12 km s and the upper limit log 
M acc < -9.67 M /y is not stringent for an object of 0.3 M , but, 
combined with the Ha properties, it suggests that this is a low 
accretor (if any). 

S0818. Classified as TD. This star shows significant excess 
emission, in the higher quartile at 3.6, 4.5, 5.8 and at 24 pm. 
We measure M flfc =-9.45 M /y, typical of optically thick disks. 
The Ha is broad (10% width of 332 + 25 km s~') and shows an 
inverse P-Cygni profile, with the emission peak at the line center 
position. The red-shifted absorption goes below the continuum 
(type IVR), confirming the evidence of a high accretion rate. 

S0897. Classified as TD; the IRAC excess emission is 
clearly detected, but lower than the <x Ori medians, while the 24 
pm excess is strong. We measure a lower limit to M acc {> -9.13 
M /y), but the star is clearly accreting. We have two measure- 
ments of the Ha profile, one with Giraffe in October 2004 and 
one with SARG acquired in January 2009. Within these two 
epochs the maxima change in position and strength, with the pri- 
mary one blue-shifted in 2004 and red-shifted in 2009, and the 
secondary one with opposite behaviour. The wavelength separa- 
tion of the blue and red emission peaks decreases from 2004 and 
2009, while the central reversal seems to be at the line center in 
the first epoch and then slightly blue-shifted in the second epoch. 
The profile changes from IIR to IIB, following the Reipurt et 
al. (1996) scheme, where these types are characterized by sec- 
ondary peaks exceeding half the strength of the primary peaks, 
as we observe. This is a common phenomenon in accreting T 
Tauri stars, probably due to the interplay of variable accretion 
and mass-loss. 

SO908. EV disk, with excess emission within the lower quar- 
tiles at all IRAC wavelengths, and significant excess at 24 pm. 
We measure M acc --931 M /y, typical of optically thick disks. 
The Ha is broad and asymmetric, with less emission in the red 
than in the blue. This type of profile (IIIR) is the less frequent in 
the scheme classification of Reipurt et al. (1996). Following the 
radiative transfer models developed by Kurosawa et al. (2006), 
this profile morphology requires some obscuration by the dusty 
disk, i.e. a high inclination, explaining the rarity of the profile. A 
highly inclined disk is consistent with the SED properties. 
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Table C.l. Accretion properties of the stars shown in Fig. IC.2I "Stars names, following the nomenclature of Hernandez et al. 2007; 
''values of the mass accretion rates derived from the U-band excess emission; c width of the Ha line at 10% of the peak; ^values of the 
M acc obtained using the empirical relation which involves the Ha width at 10% of the peak (Natta et al. 2004); c pseudo-equivalent 
width of the Ha line. 
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Fig. C.l. SEDs of four EV disks and 2 TD. The squares shows 
the observed fluxes (see Table IC.2I ). The lines plot the model 
atmosphere from Allard et al. (2000) at the appropriate T e ff, 
log g = 4.0, normalized to the J band. Each panel gives the name 
of the star (see Table IC. ft and the mass accretion rate. 
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14.87 


12.99 


12.36 


12.05 


11.77 


11.64 


11.65 


11.69 




10.962 


10.035 


9.526 


9.382 


9.2 


9.21 


9.19 


9.15 


9.03 


15.06 


13.473 


12.772 


12.528 


12.29 


12.23 


12.17 


12.25 






10.969 


10.287 


10.082 


10.05 


9.98 


9.93 


9.93 




14.18 


12.92 


12.28 


12.03 


11.82 


11.81 


11.72 


11.8 




15.12 


13.4 


12.76 


12.44 


12.07 


11.91 


11.61 


10.99 


8.58 


15.52 


13.68 


13 


12.73 


12.47 


12.5 


12.38 


12.39 






11.506 


10.879 


10.658 


10.55 


10.52 


10.47 


10.47 




14.44 


12.825 


12.064 


11.59 


10.87 


10.39 


10.13 


9.22 


6.36 


15.3 


13.46 


12.87 


12.67 


12.16 


11.86 


11.39 


10.02 


5.6 


17.184 


14.75 


14.21 


13.8 


13.51 


13.33 


13.33 


13.34 


9.8 


11.484 


10.133 


9.491 


9.293 


9.17 


9.27 


9.13 


9.12 


9.18 


12.76 


11.054 


10.251 


9.832 


9.26 


8.89 


8.56 


7.66 


5.07 




12.255 


11.568 


11.379 


11.21 


11.22 


11.22 


11.16 




13.76 


12.27 


11.51 


11.29 


11.12 


11.18 


11.03 


11 




15.17 


13.31 


12.69 


12.36 


12.1 


12.03 


11.97 


11.92 
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11.471 


10.008 


9.366 


9.189 


9.2 


9.25 


9.22 


9.01 


9.09 


S01353 


85.02959 


-2.54577 


in 


20.66 


19 


18.15 


17.9 


15.28 


13.42 


12.81 


12.54 


12.23 


12.15 


12.19 


12.1 




S01359 


85.03622 


-2.54536 


in 


19.47 




17.16 




13.9 


11.77 


11.15 


10.85 


10.59 


10.51 


10.48 


10.48 




S01361 


85.03713 


-2.55938 


11 


<17 


16.2 


15.1 




13.907 


11.501 


10.546 


9.911 


9.1 


8.72 


8.33 


7.5 




S01362 


85.03891 


-2.41853 


11 


19.8 




18.27 




15.15 


13.15 


12.5 


12.15 


11.51 


11.09 


10.66 


9.87 





